# NOT RUN {
base = iris
names(base) = c("y", "x1", "x2", "x3", "species")
# Multiple estimation
res = feols(y ~ csw(x1, x2, x3), base, split = ~species)
# The first estimation
res[[1]]
# The second one, etc
res[[2]]
# }
Run the code above in your browser using DataLab